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A striking experimental observation in vibrated granu- 
lar materials is the Brazil-nut effect (BNE), where a large 
intruder particle immersed in a vertically shaken bed of 
smaller particles rises to the top, even when it is much 
densei^ The usual practice, while describing these ex- 
periments, has been to use the dimensionless acceleration 
r = auj'^/g, where a and lu are respectively the ampU- 
tude and the angular frequency of vibration and g is the 
acceleration due to gravity Considering a vibrated 

quasi-two-dimensional bed of mustard seeds, we show 
here that the peak-to-peak velocity of shaking v ~ au, 
rather than F, is the relevant parameter in the regime 
where boundary-driven granular convection is the main 
driving mechanism. We find that the rise-time r of an 
intruder is described by the scaling law r ^ (w — Wc)^, 
where Vc is identified as the critical vibration velocity for 
the onset of convective motion of the mustard seeds. This 
scaUng form holds over a wide range of {a,uj), diameter 
and density of the intruder. 

Besides its obvious technological importance, granular 
matter is one of the most interesting examples of a driven dis- 
sipative system. Shaken granular matter exhibits various in- 
teresting phenomena such as convection, segregation and jam- 
ming. The BNE involves the motion of a single intruder in a 
vibrated bed of granular material. A related effect is the seg- 
regation obtained on vibrating a mixture of particles of differ- 
ent sizes and densities. While some progress has been made, 
a complete understanding of these observations is still lack- 
ing. The consensus now is that different segregation mech- 
anisms such as convection, void filling, air-drag and kinetic- 
energy driven methods are responsible for the effect, each be- 
ing dominant in different parameter regimes^^ In this 
Letter, we show that in a broad parameter regime, convec- 
tion is responsible for BNE in quasi-two dimensions. It is 
well-known that granular convection can be either boundary 
wall-driven or buoyancy-driven. Boundary-driven convection 
is unique to granular systems and refers to the case where the 
surface properties of the side walls and their shape are cru- 
cial in determining the formation and nature of the convection 
Buoyancy-driven granular convection, on the other 
hand, is analogous to fluid convection and is seen under strong 
driving'^' 

One of the simplest measurements in a granular convect- 
ing fluid is that of the rise-time r of a large intruder This 
quantity effectively provides a single observable characteriz- 
ing the complex velocity field in the system. The intruder rise- 
time has been used as a probe for understanding the effects of 
driving frequency and amplitude, interstitial air-flow, the size 
and density of the intruder, and wall effects on BNE. Duran 
studied BNE in a quasi-two-dimensional bed of alu- 
minium beads. They demonstrated two different segregation 




FIG. 1 : Experimental setup, a. The initial position of an acrylic 
intruder of diameter 1.7 cm embedded in a column of mustard seeds 
of sizes between 1.5 mm and 2 mm. The rise-time r is the time taken 
by the centre of the intruder to move from the height /12 to the height 
hs. A typical trajectory of the rising intruder is shown by the vertical 
white line. b. Convection pattern obtained by superposing many con- 
secutive time-lapsed images of the vibrated system. The two loops 
(yellow) with arrows indicate the flow direction, c. Intruders used in 
the experiment with relative densities ranging between 0.30 - 2.34. 



mechanisms: arching effects in the range Fc < F < 1.5 with 
Fc ~ 1, and a convective regime for F > 1.5. In the arching 
regime, intruder rise was observed only when intruder-to-bead 
diameter ratio was large enough. In the convective regime, 
the rise-time was independent of intruder diameter For a 
three-dimensional bed of glass beads. Knight et alP studied 
the rise time of tracer particles in the boundary-driven con- 
vection regime. For a fixed driving frequency they reported 
a r ^ (F — Fc)^^'^ dependence with Fc ~ 1.2, while for 
fixed F values, they obtained an exponential dependence on 
frequency. In a later work, Vanel et al} looked at the mo- 
tion of an intruder and reported deviations from the above be- 
haviour However, none of these studies have obtained any 
scaling of the rise-time data. An understanding of the depen- 
dence of rise-time on driving amplitude and frequency there- 
fore remains incomplete. 

The present work reports experimental measurements of 
the rise-time of an intruder in a vertically vibrated quasi- 
two-dimensional bed of mustard seeds (Fig. [T]i. The mea- 
surements were made under a fairly wide range of driving 
parameters, obtained by varying frequencies and amplitudes 
(2.5 ^ F < 25), intruder densities and diameters. The most 
remarkable finding of this work is the collapse of the rise- 
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FIG. 2: A series of snapshots of the vibrated bed taken at regular intervals of 8 seconds . Top and bottom panels are for experiments done 
without and witli an intruder respectively. The shaking parameters were set at / = 50 Hz, a — 1.6 mm and pr = 1.07. The alternate black 
and yellow mustard layers destabilize at the walls and eventually lead to bulk convection rolls. An inspection of these figures confirms that the 
granular convection patterns do not change significantly as a result of the presence of the intruder, except close to the intruder boundary. 



time data when plotted as a function of the driving velocity 
V = auj. Throughout the parameter regime where the col- 
lapse holds, we find that boundary driven convection is the 
dominant mechanism responsible for the rise of the intruder. 
Shaking the bed vertically results in a thin downward mov- 
ing layer of mustard seeds adjacent to each side wall, which 
sets up bulk upward motion elsewhere in the system (Fig.[T]3). 
This convective motion of the medium carries the intruder to 
the top of the bed. There is no dependence of the rise-time on 
the density of the intruder, with the intruder behaving approx- 
imately like a mass-less tracer particle tracking the convective 
flow of the background granular medium. This is also evident 
from the very similar flow patterns of the medium obtained 
with and without the intruder (Fig.[2]l. We note that v has been 
used as a scaling parameter in the vibro-fluidized regim^^^^G^l 
that is typically realized at much higher driving (F > 50) and 
is not wall-driven. However, no scaling has been achieved so 
far in the boundary-driven convection regime. Given the non- 
homogeneous flow fields that are set up in our experiments, 
the dependence of r on shaking velocity is rather surprising. 

In our experiments, we measure the rise-times r 
for several peak-to-peak amplitudes a and frequencies 
/ {— Lo /2tt), for various relative densities = 
(intruder density) /(mustard seed density). The rise-time 
decreases with increasing shaking frequencies and ampli- 
tudes, as seen in the upper inset of Fig. |3] It is observed 
that r is almost independent of density over the range used 
in these experiments (pr = 0. 30 - 2 .34). This is in contrast 
to the three-dimensional result J^E^Ell where a non-monotonic 
dependence of r on pr was found. There is no scaling of r 
with r (lower inset of Fig. |3]l. In contrast, a plot of r as a 
function of the shaking velocity v = auj (Fig.[3]l shows an ex- 
cellent collapse of the data. This collapse holds for the entire 
range of intiTider densities and the two intruder sizes (diameter 
1.7 cm and 1 cm) that were used. The collapsed r-vs-u data of 
Fig.plis well described by the scaling-law t — A{v ~ Vc)^'^, 




20 30 40 50 60 70 

V = flco (cm/sec) 



FIG. 3: Intruder rise-time plots. Upper inset: Plots of the rise-time 
T versus shaking amplitude a for different frequencies: i) 20 Hz (□) 
ii) 35 Hz (A) and iii) 50 Hz (Q) for an intruder of diameter 1.7 cm 
and for four different relative densities pr = 0.3 (blue), 0.8 (green), 
1.07 (red) and 2.34 (magenta). The filled symbols correspond to a 
1 cm diameter intruder with pr — 1.07. Lower inset: Plot of the 
same data versus the dimensionless acceleration F. Main figure: 
Plot of r versus peak-to-peak velocity of shaking v — auj. The data 
collapses for the entire range of shaking velocities, for all relative 
densities and for the two intruder sizes. The solid line shows the 
functional form r ~ (w — Wc)"^ with = 16 cm/sec. 



where A and Vc are fitting parameters. For the range of pr 
values investigated, Vc is found to be around 15-18 cm/sec 
and A is around 13, 000 ± 500 cm^/sec. 

The absence of any dependence of the rise-time on the size 
and density of the intruder implies that convection is the dom- 
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FIG. 4: Convection velocities at tlie side walls. The convection 
velocities, in the absence of the intruder, are plotted versus the peak- 
to-peak velocity of shaking v = auj for shaking frequencies: i) 20 
Hz (□), ii) 35 Hz ( ) and iii) 50 Hz (Q). Each linear fit (denoted 
by solid line) is extrapolated to estimate the peak-to-peak shaking 
velocity for the onset of convection. The x-intercepts lie in the range 
of 16 - 17 cm/sec. Inset: Plot of the convection velocities at the 
side-wall as a function of F. 



inant driving mechanism. To confirm this premise, we make 
an independent measurement of the convection velocity of the 
mustard seeds without the intruder. Over a range of shak- 
ing frequencies and amplitudes, we measure the velocity of 
tagged mustard particles which are moving down along the 
two side walls. The velocity of the boundary layer particles is 
an order of magnitude larger than the mean velocity of the in- 
truder This, when plotted against v = aoj, shows a much bet- 
ter data-collapse than when plotted as a function of F (Fig.[4|. 
A linear extrapolation of the data to zero convection velocity 
gives a critical velocity of 16 - 17 cm/sec for the onset of con- 
vection. This critical velocity range overlaps with the range 
of the fitting parameter Vc obtained from the rise-time data. It 
is therefore natural to identify Vc as the critical shaking veloc- 
ity required for the onset of convection cycles in the granular 
column. 

Our understanding of the different parameter regimes for 
boundary driven convection is summarized in the phase- 
diagram in Fig. [5] We note that while F > 1 is necessary 
for bulk convection, it is not a sufficient condition. As seen 
in the phase diagram, there is a region of parameter values, 
with F > 1 and v < Vc, where no bulk convection occurs. 
In this regime, we observe some motion of the mustard seeds 
restricted entirely to the the upper corners of the bed. Such 
conditions cannot give rise to BNE for the intruder sizes stud- 
ied here. We expect that for much larger intruder sizes, other 
mechanisms may come into play which can drive BNE^. We 
also find that on varying the bed-height, the scaling of rise- 
time continues to hold for smaller heights (supplementary ma- 
terial) but breaks down for larger hegihts. The parameters 
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FIG. 5: Phase diagram. Each data point (o) corresponds to a set 
of (amplitude, frequency) parameter values used in our experiments. 
The solid line corresponds to w = a{2nf) = Vc, while the dashed 
line corresponds to F = a(27r/)^/(/ — 1. Bulk convection occurs 
only in the region where both the conditions F > 1 and v > Vc are 
satisfied. 



A and Vc depend on the bed height. The microscopic ori- 
gin of the parameters A and v^, is not clear. They depend on 
the bed-particle diameter and the bed-height, which are eas- 
ily characterizable quantities, but also on details such as the 
wall, and bed-particle roughness, which are harder to char- 
acteize. Finally, we anticipate that a similar scaling should 
hold even in three-dimensional granular systems, in some pa- 
rameter regime. 

Methods: The experimental setup consists of a quasi-two- 
dimensional rectangular acrylic cell of height 25 cm, width 13 
cm and gap 6 mm. The cell is filled up to a height h with black 
mustard seeds of density 1.11 g/cc and sizes between 1.5 and 
2 mm. The intruder is a large acrylic disk of diameter 1.7 cm, 
thickness 5 mm and density 1.07 g/cc. The thickness of the 
intruder is chosen such that mustard seeds are not able to enter 
the space between the intruder and the cell walls. The cell is 
mounted on a 50 kgf DESPL electrodynamic shaker executing 
vertical sinusoidal vibrations of the required amplitudes and 
frequencies. The narrow side walls of the cell are roughened 
by lining them with fine sandpaper (ISO grit designation P — 
100). For these boundary conditions, the experiments were 
reproducible. However, for smooth walls it was found that the 
convection patterns were irreproducible and often the system 
would jam. The humidity of the lab was maintained at 40%. 

In order to produce fixed initial conditions corresponding 
roughly to a fixed spatial distribution of mustard seeds, the 
cell was initially vibrated at a high F and the intruder was 
pushed to a fixed depth hi (Fig. [T^). The driving was then 
removed and the top layer leveled. The cell was then vi- 
brated at the required a and /. The upward motion of the 
intruder was recorded using a Panasonic Lumix DMC-TZ3 
digital camera at a rate of 33 frames/sec. High speed track- 
ing of mustard-seed convection at the side walls was done 
using an IDT MotionPro high speed camera at a rate of 150 
frames/sec. The tagged particles for this experiment were yel- 
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low mustard seeds of sizes and densities within the same range 
as the black mustard seeds. These were dispersed randomly 
in the bed. The tracking of the intruder (Fig. [T^) and the 
coloured particles near the side walls was done using a Lab- 
view based particle tracking algorithm^". In order to change 
the density of the intruder, several acrylic intruders with cen- 
tral holes were made. By changing the size of the air holes or 
inlaying them with copper, a range of density ratios p,. = 0.30 
- 2.34 (Fig.[T]:) was achieved: pr — 0.3 for acrylic with 1.4 
cm air hole, 0.8 for acrylic with 1.2 cm air hole, 1.07 for solid 
acrylic disk and 2.34 for acrylic with 0.65 cm copper 

The rise-time of the particle was computed by estimating 



the time taken by the intruder to rise from a fixed starting po- 
sition h2 to a fixed final position h^^ (Fig.[T^). The mean rise- 
time T was computed by averaging over 10 trials for every 
parameter set. For the convection velocity at the wall, the fall 
time of only those seeds which remained in the first boundary 
layer were considered and the mean velocity was obtained by 
averaging over many seeds. 

The granular convection phenomenon was studied in 
greater detail by filling up the rectangular cell with alternat- 
ing layers of black and yellow mustard seeds, with each layer 
being approximately 2 cm thick. 
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FIG. 6; Intruder position vs. time plot. Position vs. time plots for an intruder of pr — 1.07, for i) / — 20 Hz, a = 4.8 mm (green), ii) 
/ = 35 Hz, a = 2.0 mm (red) and / = 50 Hz, a — 1.2 mm (blue). The dashed horizontal lines (black) represent /i2 and hs (see Fig. la of 
main paper). The rise-time of the particle is computed by estimating the time it takes for the intruder to rise from the position /12 to /13. The 
vertical dashed lines in green, red and blue enclose the time windows required by the intruder to rise between h2 and /13 for the three shaking 
experiments (i), (ii) and (iii) respectively. 
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FIG. 7: Different realizations of intruder trajectories. Position vs. time plots for 10 different trials of the rise of a p, = 1.07 intruder for 
/ = 50Hz and a = 1.4 mm. The trajectories show roughly quadratic dependence on time and different repeats of the experiment for fixed 
control parameters lead to slightly different trajectories and correspondingly different rise- times. 
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FIG. 8: Plots of the rise-time for different bed-lieiglits. Plot of r versus v = auj for two different bed heights h = 8cm (green) and h — IQ 
cm (red), for three different frequencies 20 Hz (□), 35 Hz (A) and 50 Hz (O)- In all experiments, the relative density was p,. = 1.07 and 
the intruder size was 1.7 cm. We see good collapse of data for both heights. The fits are to the scaling form A/{v — Vc)^ with A ~ 13000 
cm^/sec, Vc ~ 15 cm/sec for the 10 cm bed and A ~ 7000 cm^/sec, iic ~ 7 cm/sec for the 8 cm bed 
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(a) / = 20 Hz, a = 3.7 mm, auj = 46.5 cm/sec 




(b) / = 35 Hz, a = 2.1 mm, auj = 46.2 cm/sec 




(c) / = 50 Hz, a = 1.5 mm, auj = 47.1 cm/sec 

FIG. 9: Tagged particle trajectories (y-vs-t) near wall. Each bright, slanted yellow trajectory represents the time-evolution of the y (vertical- 
coordinate of a tagged particle (yellow mustard seed) adjacent to the side wall, for three different parameter sets. The slope of the yellow 
trajectories give the velocities of the tagged particles as (a) 4.2 cm/sec, (b) 5.0 cm/sec and (c) 4.3 cm/sec respectively. The white sinusoidal 
trace (horizontal) represents the time-evolution of a reference white mark on the acrylic cell. 



